Impurity in a Luttinger liquid away from half-filling: a numerical study 
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Conformal field theory gives quite detailed predictions for 
the low energy spectrum and scaling exponents of a massless 
Luttinger liquid at generic filling in the presence of an impu- 
rity. While these predictions were verified for half-filled sys- 
tems, there was till now no analysis away from this particular 
filling. Here, we fill in this gap by numerically investigating a 
quarter-filled system using the density matrix renormalization 
group technique. Our results confirm conformal field theory 
predictions, and suggest that they are indeed valid for arbi- 
trary fillings. 

PACS numbers: 71.10.Pm, 72.15Nj 



For a half-filled Luttinger liquid, these predictions have 
been numerically verified (^^]. In the present work, we 
study numerically the non-half-filled system, which is the 
generic situation for realistic systems. Compared with 
the half-filled case, the non-half-fillcd one has non-zero 
forward scattering phase shift at the Fermi surface, which 
makes the analysis more cumbersome. In spite of that, we 
show that, also away from half-filling, the correct fixed 
point boundary conditions are those predicted in Refs. 




At the end we will discuss in detail the boundary con- 
ditions in relation to some previous studies || which gave 
results different from those mentioned above. 



I. INTRODUCTION 

The behavior of a one-dimensional (ID) interacting 
Fermi system (Luttinger liquid) in the presence of a sin- 
gle impurity has been the subject of an intensive the- 
oretical and numerical investigation in recent years, for 
its interesting anomalies with respect to higher dimen- 
sions, and its implications to a variety of physical prob- 
lems as for instance the behavior of quantum wires (see 
e.g. Refs. [|l|,||) or the tunneling through a constriction 
in the fractional quantum Hall regime ||. Specifically, 
if the interaction is repulsive, the electrons at low ener- 
gies feel the impurity potential as if it were effectively 
infinite. On the contrary, for attractive interactions, the 
effective scattering close to the Fermi energy is vanish- 
ingly small (see e.g. Ref. [[|). In their seminal paper, 
Kane and Fisher [Q argued that the low energy behavior 
for repulsive interaction in fact corresponds to a chain 
disconnected at the impurity site. In the renormaliza- 
tion group language of boundary conformal field theory 
(CFT) (see e.g. Refs. H|j), f° r repulsive interaction the 
open boundary condition (OBC) is the stable fixed point 
while the periodic boundary condition (PBC) is unstable, 
and oppositely for attractive interaction (except for the 
presence of nontrivial fixed points for spinning electrons 
with spin anisotropic interaction [7p^]). CFT provides 
the low energy spectrum and the scaling exponents both 
for the approach to the appropriate stable fixed point and 
for the departure from the unstable one. 



II. NUMERICAL CALCULATION 

We study numerically the S — 1/2 XXZ Heisenberg 
chains with anisotropy a = J z / J xy = ±0.5 at magnetiza- 
tion per site M — 1/4, which correspond to quarter-filled 
Luttinger liquids with repulsive (a = 0.5) or attractive 
interaction (a = —0.5). We model the impurity by mod- 
ifying the strength of a particular bond, scaling it by a 
factor b. Specifically, we consider the following Hamilto- 
nian 

L-l 

H = [(S^S i+1 + S t ± )/2 + aSf Sf +1 ] p 

+ 6 [(S+S£ + S^S+)/2 + aSfSi] , 

where and S~ are spin raising and lowering opera- 
tors at site i, Sf being its z-component. L is the chain 
length. 6 = and b — 1 correspond to what we will de- 
fine as ideal OBC and PBC , respectively. We use the 
density matrix renormalization group (DMRG) method 
ifiof to calculate the low energy spectrum and overlaps 
between ground state wave functions. To ensure the cor- 
rect magnetization M = 1/4, we insert four sites at each 
step of DMRG iteration. We run DMRG for a = ±0.5, 
L = 8, 12, . . . , 60, and b = 0.0, 0.05, 0.1, . . . , 0.9, 0.95, 1.0. 
The optimal number of states which are kept is 300 and 
the truncation errors are less than 5 x 10 -7 . Like our 
previous DMRG calculation for magnetization M = 



1 



||, we calculate the lowest energy states with parity 
even and odd (with respect to the modified bond), and 
z component of total spin 5'°' = ML and ML + 1 for 
M = 1/4. We denote these energies in increasing order 
as Eq, Ei, E2, and E3. For states with magnetization M, 
Eq corresponds to the ground state energy while Ei cor- 
responds to the energy of the lowest excited state with 
different parity. E\ and E3 correspond to states with 
S l z ot = ML + 1, but with different parity. In addition 
to these low energy levels, we calculate also the ground 
state wave function overlaps between PBC chains (6 = 1) 
and impurity chains (b =/= 1). 



III. ENERGY SPECTRUM AND SCALING 
PREDICTED IN CFT 

The spin chain at non-zero M and with a modified 
bond can be written in terms of a boson field, 0, which is 
defined periodically: = + 2irR, and of its conjugate 
momentum II. A detailed derivation is given in Ref. §] 
for the case M = 0. Here we just point out the differences 
that occur for non-zero M. The compactification radius, 
R, and the spin- wave velocity, v s , depend both on a and 
on M. Their values for a = ±0.5 and M = 1/4 are given 
in Table I. The effective Hamiltonian is: 
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(3.1) 



[We will generally set the spin-wave velocity v s to 1 in 
most of what follows.] Here the modified bond is parame- 
terized in terms of two boundary operators, with coupling 
constants Vf and Vf,. The Vf term is allowed, and there- 
fore expected to occur, whenever S z — > — S z symmetry 
as well as translational invariance are broken, which is 
the case for M ^ and with a modified bond. Specif- 
ically, Vf and Vb correspond to forward and backward 
scattering in the Luttinger liquid, respectively. Forward 
scattering is marginal but backward scattering is relevant 
(irrelevant) for R > l/VSvr (R < I/v^tt). 

In cases where backscattering is irrelevant, which we 
define generically as PBC cases (not to be confused with 
ideal PBC with b = 1), Vf induces a discontinuity in 
at the impurity site, 0(0+) - 0(0-) = 0(0) - 0(L) = 
2R5 modulus 2txR. Therefore the mode expansion can 
be written as: 
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(3.2) 



Here the dots represent the harmonic oscillator modes, 
ri = S^ ^ — ML and m arc integer quantum numbers. 
(The periodicity is preserved in time by the discrete 
quantization of the momentum conjugate.) Thus the 
finite-size spectrum is: 



E = cqL + e\ + n/z 
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(3.3) 



ra^ ,R are occupation numbers for bosonic excitations of 
momentum ±2im/L. The chemical potential, fi, is re- 
lated to the external magnetic field necessary to induce 
the magnetization, M, in the spin chain. S is the anal- 
ogous to a forward scattering phase shift. Here we have 
included the zero point energy — 2nv s /(12L) which was 
ignored in Ref. Jo] . This formula differs slightly from 
the one in Ref. |q] because in that case spinless fermions 
were considered rather than S — 1/2 spins. While the 
models are identical under the Jordan- Wigner transfor- 
mation, PBC on the spins do not correspond to PBC on 
the fermions. 

When the backscattering term is relevant, 0(0) gets 
pinned. In analogy with the previous case, we define this 
situation as OBC, again not to be confused with the ideal 
OBC (6 = 0). Without forward scattering (for M = 0) 
0(0) will simply be pinned at a unique value, 0(0) = 
assuming Vb > 0. With both backward and forward scat- 
tering present, we expect to be pinned at equal and 
opposite values, ±0o at x — * O*, with a discontinuity 
at x = 0. Pinning forces m = in the mode expan- 
sion, changes the wavef unctions of the harmonic modes 
from e ±l27Tnx / L to sin(7rna;/L) and reduces the zero point 
energy by a factor of 1/4. Thus 



E = e L + ei + rifi H — — 
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The parameters eo, V s , R and /1 being bulk properties 
are the same for PBC and OBC cases. The parameters 
ex and 6, on the other hand, depend on the boundary 
conditions, as well as the 1/12 and 1/24 terms. The 
latter can be calculated by elementary methods for the 
XX spin chain with PBC or OBC, respectively. They 
have also been calculated, using the Bethe ansatz for the 
general XXZ model. They follow more generally from 
conformal invariance. El 12 Note that forward scatter- 
ing always introduces a discontinuity in at the origin, 
with or without pinning of depending on whether or 
not backscattering is relevant. However, we expect the 
value of the discontinuity cx S to depend on whether or 
not pinning occurs. 

The effect of a parity transformation on the various 
eigenstates were worked out in Ref. In the periodic 
case it takes m — > — m, <-> m R and multiplies wave- 
functions by (— l) n+m . In the case of OBC it multiples 
wavefunctions by (—1)™ and transform a bosonic opera- 
tor a m , corresponding to an excitation with wave func- 
tion sin(7rma;/L), into (— l) m a m . 
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We are now in position to work out the quantum num- 
bers and energies of the four lowest energy states. Con- 
sider first the states with Sl ot = ML. These states have 
ri = 0. The ground state has all other quantum numbers 
as well. For PBC, a state of reversed parity can be 
obtained cither by t akin g m = 1 (and all other quantum 
numbers 0) in Eq. (3.3) or by taking an odd parity lin- 
ear combination of the states with (mf,mf) = (1,0) 
or (0,1) (and all other quantum numbers 0). These 
two states have energies measured from the ground state 
energy of (2ttv s / L)/ (AttR 2 ) and (2ttv s /L) respectively. 
Thus which of these states has lower energy depends on 
whether AttR 2 > 1 or AttR 2 < 1. We see from Table I 
that the first one has lower energy for a = 1/2, while the 
second one has lower energy for a = —1/2. The states 
with Sl ot = ML + 1 have n = 1. The state with all other 
quantum numbers has reversed parity relative to the 
grounds state. For PBC the state with the same parity 
as the ground state is obtained again by choosing m = 1 



the correction to d 2 — do becomes larger than the fixed 
point value, at sufficiently large L. Since the state |2 > 
has m = 1, and m corresponds to a conserved quantum 
number at the PBC fixed point, the only non-zero matrix 
element of cos(cj)/R) in the subspace of |0 > and |2 > is an 
off-diagonal one < 2| cos(<fi/R)\0 >. Since this is larger 
than the diagonal matrix elements, the energy difference, 
divided by L, scales as L v , where y is the renormaliza- 
tion group eigenvalue of the associated coupling constant, 
y = 1 - x. On the other hand, for AttR 2 < 1, the PBC 
fixed point is stable, and in fact the corrections to the 
di's from O vanish at large L. In this case, the state 
1 2 > has m = (and m\ or — 1). Hence the first 
order corrections to the energies vanish and the leading 
corrections are second order, that is d 2 — do oc 

At the OBC fixed point, the most relevant permit- 
ted operator is a product of the operators which change 
ri ± 1 from the two sides of the broken chain, that is 



S 1 S L 



The n = ±1 operator has dimension 



(for AttR 2 > 1) or either or mf = 1 (for AttR 2 < 1). 2ttR 2 (1 =f 28 /tt). Since the two sides of the broken chain 



For OBC, the lowest energy states for given S 1 ' ' (i.e. n) 
are obtained by choosing all the m n = 0, while the par- 
ity reversed lowest energy states are obtained by choosing 
mi = 1- We define new quantities, di, by writing the 1/L 
term in the energy as 
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Ei = ... H \-di 

L \ 12 4 



(3.5) 



for PBC or OBC. The values of the di for the first four 
states for ideal PBC and OBC for a = ±0.5 are given in 
Table II. 

The renormalization group flow in the vicinity of the 
PBC and OBC fixed points are determined by the per- 
mitted boundary operator, O, of lowest scaling dimen- 
sion, x. Depending on whether this dimension obeys 
x < 1 or x > 1, the fixed point will be unstable or stable. 
We may calculate the (size dependent) corrections to the 
energies, in the vicinity of the fixed point, by doing low- 
est order perturbation theory in this operator. |13| The 
dimensions of a boundary operator at the PBC and OBC 
fixed points can be read off from the finite size energies. 
Specifically, if a boundary operator modifies some quan- 
tum numbers, its dimension is simply given by the term 
inside the square brackets of Eq. (3.3) for PBC and (3.4) 
for OBC with the appropriate change of quantum num- 
bers, minus the analogous term with the ground state 
quantum numbers. 

For instance, the leading boundary operator which can 
be defined at the PBC fixed point is the backsc atte r- 
ing potential O = cos(4>/R) [see the Hamiltonian (3.1)1. 
Since it changes m by ±1, its scaling dimension is simply 
x = 1/AttR 2 . The lowest order perturbative calculation 
is quite different depending on whether AttR 2 > 1 or < 1. 
Let us denote |0 > and |2 > the states corresponding to 
the energy levels Eq and E 2 pre viously introduced, and 
use do and d 2 defined in Eq.(3.5). If AttR 2 > 1, then O is 
relevant and the PBC fixed point is unstable. In this case 



are decoupled at the ideal OBC fixed point, the dimen- 
sion of the product of operators on the two sides is simply 
additive, giving a total dimension for the leading opera- 
tor of AttR 2 . Since there is no extra conserved quantum 
number at the OBC fixed point we expect this operator 
to have a non-vanishing expectation value in the states 
|0 > and 1 2 >. Hence the correction to d 2 — do should 
scale as L v . 

Exactly at b = or 1, the set of allowed operators is 
reduced. For the PBC case, b = 1, all local as well as 
777 7^ operators are forbidden by translation invariance. 
The leading irrelevant operators correspond to 



dx{dcj)/dx){d(j)/dx) 2 



dx{dcj>/dx){dcj>/dt)' 1 



(3.6) 



of dimension 2. [d<j)/dx can be eliminated by a redefi- 
nition of (f).] These dimension 2 operators are forbidden 
by symmetry for M = and so did not appear in [^). 
They contribute 1/L 2 corrections to the energies. For 
the 6 = case, since translational symmetry is broken, 
we can construct boundary operators but only at one side 
or the other of the broken link; we can not take a product 
of operators on both sides as above. In addition, spin 
rotational invariance forbids all operators with 77 ^ 0. 
Again d(j)/dx is eliminated by a redefinition of 4>, so the 
leading irrelevant boundary operator is (d<fi(0) / dx) 2 of 
dimension 2 This also contributes 1/L 2 corrections 
to the energies. 

When b is close to or 1, we obtain analytic 1/L 2 
finite-size corrections, as just discussed, with coefficients 
of 0(1), as well as non-analytic finite size corrections, 
as previously shown, with coefficients that vanish in the 
limit b — > or 1. 
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IV. NUMERICAL STUDIES ON SPECTRUM 
AND RENORMALIZATION 

We first consi der t he bulk quantities, eo, fJ>, v s and R, 
denned in Eq. (3.3), We determined these using both 
the DMRG and the Bethe ansatz technique |]l6| . 

These bulk quantities can be determined from the 
finite-size spectrum with PBC which can be derived 
from the Bethe ansatz |Lffl. The structure of the spec- 
trum up to 0(1/ L) term agrees with the CFT prediction 



[S = 0, e% = in Eq. (|3.3|) ]. By comparison, the various 
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Here the right-hand sides are given by 
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K(rj - rf)p(rf)drf 



: -o(v)p(v)dv, 



where cos0 = a (see Eq. (2^1)) and the integral kernel is 
given by 



K(rj) = 



tan0 



tan 2 cosli (rj/2) + sinli (77/2) 



(4.3) 



The parameter A is fixed from the filling factor v = 1/4 
by 



a(r])d-q = v, 



(4.4) 



and the chemical potential fi is determined by the condi- 
tion e,j(A) = 0, where 

td(ri) = eo(v) ~ i: I K(rj - V ')€ d (r,')dr,' . (4.5) 



2tt 



-A 



The half- filled case v = 1/2 corresponds to A = 00. In 
this case the integral equations can be solved analytically 



by Wiener-Hopf method jl7) , which gives R = y 2 J n 
and v s = 7r sin 0/(29), as deduced earlier from the Bethe 



ansatz solution using other methods. [ |18|Jl9| In partic- 
ular, for a = ±0.5 and M — 0, the case considered in 
Ref. H, v s = 3V3/4, 3^3/8, and AirR 2 = 2/3, 1/3, 
respectively. The Bethe ansatz results for L(E\ — Eq), 
e 2 = (E 2 - E )/(E 1 - E ) and e 3 = (E 3 - E )/(E 1 - E ) 
are given in Table III. Excellent agreement is obtained 
with the finite-size spectrum calculated using DMRG in 

For the non-half-filled case, these equations cannot be 
solved analytically. We solved them numerically by iter- 
ative numerical integration. Since we have already taken 
the limit of large system size L in analytic treatment, we 
can obtain highly accurate results with relatively small 
amount of computation. Accuracy to six digits can be 
achieved within minutes on an IBM RS6000 workstation. 
On the other hand, it would be difficult to calculate other 
quantities, for example orthogonality exponent, directly 
from Bethe Ansatz. 

By fitting all f our Ej 's o btained from DMRG accord- 
ing to Eqs. ( |3.3[ ) and (3.4), for the cases we defined as 
ideal PBC (b = 1) and OBC (6 = 0), we find consis- 
tent bulk quantities, eo, /x, v s , and R, which are listed in 
Table I. These values agree well with the Bethe Ansatz 
results, as also shown in Table I. Thus the accuracy of 
our DMRG calculation is confirmed. We also determined, 
from DMRG, the non-zero forward scattering phase shift, 
6, for the case of OBC, as shown in Table I. 

With these parameters, we have found that the finite- 
size corrections to the energies, di, found numerically 
agree with the low energy spectrum predicted by CFT 
for PBC and OBC and listed in Table II. The finite size 
values of di are plotted versus 1/L in Figures 1 and 2, 
and extrapolated at L — > 00 by fitting with polynomials 
of 1/L . The quality of these fittings confirms the pre- 
diction that only analytic finite size corrections exist in 
the cases with b = or 1. 

After having verified the low energy CFT spectrum, 
we confirm another prediction, namely that the orthogo- 
nality exponent between the ground states of two models 
described by the same Hamiltonian but different confor- 
mally invariant boundary conditions, is equal to the dif- 
ference between their corresponding cio's. The orthogo- 
nality exponent a is defined through the overlap between 
the two ground state wave functions, \<p) and \4>o), for a 
chain of size L. In particular, 



0(L) = 



(4.6) 



The overlap integrals O(L) are evaluated numerically. 
The a for L — > 00 limit is obtained by extrapolating 



a(L) = 



lnO(L + 4) -hiO(L) 
lni-ln(L + 4) 



(4.7) 



We start by verifying the above prediction for the overlap 
between the ground states with ideai PBC and OBC. 
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The ideal OBC d 's for L -> oo limit (for the idea] PBC, 
c?o = 0) is obtained by extrapolating 



d (L) = 



1 L(L + 4) 



12 



87TV, 



[ J B (i + 4)-£o(i)-4eo]. 

(4.8) 



In Fig. 3 we plot a(L) 7 the exponent of the overlap be- 
tween the ideal OBC and PBC ground state wavefunc- 
tions, and the ideal OBC do(L), vs. l/L. Up to the two 
significant digits the extrapolated a and d a arc equal for 
L — ► oo limit, for both a = 0.5 and a = —0.5. 

We point out here that the low energy spectrum in 
Table II holds also if the magnetization M — 0, as it was 
shown in Ref. [^J. For generic impurity strength, it was 
shown in that reference that the M = system flows 
either to OBC (a > 0) or to PBC (a < 0), by a detailed 
analysis of the size dependence of the low energy levels 
and by calculating the orthogonality exponent. Since the 
non-universal forward scattering phase shift 6 = for 
M = 0, universality is recovered and all systems flowing 
to the OBC fixed point had an extrapolated d = 1/16, 
in agreement with CFT prediction. Moreover, 1/16 was 
also shown to be equal to the orthogonality exponent of 
the overlap between the ground states with PBC in the 
presence and in the absence of the impurity. 

Analogously, in the following, we will first show that 
the low energy spectrum for M =1/4 systems does flow 
to the PBC or OBC fixed points, depending on the sign 
of the interaction. Then we will show the qualitative 
agreement between a(L) and do(L) for arbitrary impu- 
rity strength b at finite size. 

We start by showing how the low energy spectrum 
flows away from PBC towards the OBC fixed point for 
a = 0.5 and vice versa for a = —0.5. We previously said 
that 5 is not a universal property. This implies that, if a 
system flows to one of the two fixed points, it does not 
mean necessarily that the value of 5 will be exactly the 
same as that of a chain with b = or b = 1 (as given in 
Table I). Thus we expect 5 to be some non- universal func- 
tion of b and a. The finite size spectrum should be given 
by the OBC formula of Eq. and Table II for a = 0.5 
and all b but with some unknown function 5(b). Likewise, 
for a — —0.5 th e sp ectrum should be given by the PBC 
formula of Eq. ( |3.3| ) and Table II for all b but with some 
other unknown function 6(b). 8(b) should approach zero 
as b — > 1 for a = —0.5 and 6(b) should approach its OBC 
value (i=a — 0.297r) as b — > for a = 0.5 since in these two 
cases 1 — 6 or b produce no relevant operators and the 
marginal operator has a coupling constant which goes to 
zero. 

For the four energy levels we decided to calculate nu- 
merically, E 3 — Ei = E 2 — Eq always holds, as seen from 
Table II. Notice that the forward scattering phase shift 
S generated by the impurity contributes to E\ — Eq (and 
also E3 — Ei). In order to simplify the analysis, we will 
concentrate only on E2 — Eq (i.e. d 2 — do), which instead 
depends on the universal properties identifying each fixed 



point. Let us denote A(a,b) = d^a, b) — do(a,b). Then, 
for small deviations of b away from one of the two fixed 
point values 6* = 1,0, the following scaling should hold 



A(a, b) - A(a, K) cx L^°" b '\ 



(4.9) 



These exponents can be obtained by calculating the 
first non- vanishing correction to A (a, 6), within a per- 
turbation expansion in b — 6», following the prescrip- 
tions given in the previous section. The predicted 7's 
are listed in Table IV. In Fig. 4, for a — ±0.5, we plot 
ln(A(a, b) — A(a, ?>*))/ lnL vs. 1/lnL. The extrapolated 
scaling exponents 7 in the figure agree with those listed in 
Table IV, up to the two significant digits. These results 
have a very simple interpretation. If 7 < 0, the devia- 
tion from the fixed point is not a relevant perturbation 
and the system flows back to the 6»-fixed point, which is 
therefore stable. On the contrary, if 7 > 0, the perturba- 
tive correction blows up as L — ► 00, which implies that 
the system flows away from the 6* fixed point, which is 
therefore unstable. Again, our results demonstrate that, 
for repulsive interactions, the stable fixed point bound- 
ary condition is the OBC, and vice versa for attractive 
interaction. 

We have also calculated A(a, b) for arbitrary deviations 
6-6*, namely for b = 0.05, 0.15, 0.85, 0.95. To an- 
alyze the data, we made a very simple scaling ansatz in 
order just to reproduce the correct scaling behavior close 
to the OBC fixed point. Specifically, A(a, b) is plotted 
in Fig. 5 vs. arctan [L~ - 2 tan.(6?-)] /it for a = 0.5, and 
in Fig.6 vs. arctan [L - 28 tan(6f )] /tt for a = -0.5. The 
arrows in the figures show towards which fixed point the 
system flows as the chain length L increases. We see 
qualitatively how the excitation energy flows between 
OBC and PBC fixed points. The proposed scaling is 
tan [(A — |)7r] w i 7 tan(6^). Surprisingly, our scaling 
ansatz, which is in principle valid only close to b = and 
L 1 being small, seems to give a good description of the 
data for a wide range of b's. For |L 7 — ► 0, we have the 
scaling relation: 



-L 1 . 



(4.10) 



Finally we compare the orthogonality exponent a(L) 
of the overlap between the ground state with PBC and in 
the presence of the impurity, and the ground state with 
PBC but in the absence of the impurity. According to 
CFT, this exponent should also be equal to do for chain 
with impurity. We plot a(L) and do(L) vs. the same 
variables as used above for A. The resulting data are 
plotted in Fig. 7 for a — 0.5, and in Fig. 8 for a = —0.5. 
They show that a(L) and do(L) are of similar magnitude 
and have similar flows for various chain lengths and impu- 
rity strengths. Again, an important remark is needed re- 
garding the forward scattering phase shift. Since we can 
not exclude its generation at the fixed point, this would 
imply, for instance, that the orthogonality exponent for 
attractive interaction is, in general, non-zero. With long 
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enough chains it should be possible to determine 5(b) and 
check that both a(L) and do(L) extrapolate to the same 
numbers, determined by these phase shifts. However, 
with the limited lengths available (L < 60) we have not 
found it possible to do this. This is presumably due to the 
unknown value of the phase shift, the slow convergence 
of the finite-size corrections due to the small exponent 
in the non-analytic corrections (L~ 0,2 for a — 0.5) and 
the relatively large coefficient in the analytic corrections. 
The situation appeared somewhat better in the M = 
case H. 



V. DISCUSSION ON BOUNDARY CONDITION 

With periodic boundary conditions, the boson may be 
written in terms of decoupled left and right-moving parts: 



4> L (t, -x) = -<j>n(t, x) + 4> - 



(5.5) 



(j>(t, x) = 4> L (t + x) + 4> R (t - x). 



(5.1) 



The dual field is defined as 4> = 4>l — 4>R- The bound- 
ary condition, fii^) — ±0o couples left movers to right 
movers, separately on both sides of the origin: 



*k(«,0 + ) 



-0k(*,O-) 

-<M«,o+)- 



(5.2) 



On the other hand the fields on opposite sides of the 
origin are not coupled by the boundary condition. This 
means that an incoming excitation at x < is reflected, 
with unit probability with a phase shift proportional to 
(f>o and similarly for an incoming excitation at x > 0. 
This boundary condition is the same as the one used in 
Ref. S where a different basis was used. 

This coupling of left and right movers by relevant 
backscattering changes Green's functions. To illus- 
trate this, we consider the M = case where there is 
no forward scattering term, 0o = 0. One term in the 
bosonized representation of S~ is: 



i2nR4>+i(t>/R _|_ 
x e i(2-xR+l/R)<f> Le i(-2irR+l/H)4> R + 



(5.3) 



The Green's function, with PBC factorizes into separate 
Green's functions for left and right movers, giving: 

(S+(t, X)S~(0, X 1 )) CX (t + X - x ')-(2nR+l/Rr/4,n 

x(t-x + x ')-(-^R+i/R) 2 /^ + . . . 

(5.4) 

On the other hand, with the perfectly reflecting bound- 
ary condition, this correlation function equals if a; and 
x' are on opposite sides of the origin. To calculate the 
Green's function when they are on the same side (say 
x, x' > 0) we use a standard device from boundary CFT. 
p4[ We can regard the second boundary condition of Eq. 
fl5.2| ) as defining cf>R for all x > as the analytic contin- 
uation of (f>L to the negative x-axis: 



The reason this works is because the boundary condition 
is true at all t and 4>l and <j)R. on ly depend on the variables 
t + x and t — x, respectively. We emphasize that this 
analytically continued 4>l is not the original 4>l at x < 
which is completely decorrelated from the fields at x > 0. 
Using this approach, S~ (t, x) becomes bilocal: 

S~(t,x) CX e l ( 2 -xR+l/R.)<t> L (t+x)~i{-2TTR+l/R)4> L (t~x) + 

(5.6) 

Instead of getting a product of 2-point Green's functions 
for left and right moving fields, we now obtain a 4-point 
Green's function for left-movers only, giving: 



(S + (t,x)S-{0,x')) oc 



t 2 -(x+x') 2 



l/4irR 2 +TrR 2 
l/47r_R 2 +7ri? 2 



t 2 -(x-x') 2 

t — x+x 
t+x — x' 



(5.7) 



+ 



In the limit, x,x' — > + , 

S-(t,0) oc e ^fl0(*,o+) K e i4 W R0 £ (t,o+) 5 



giving: 



< S + {t,0+)S-{0,0 + ) >oc \t\ 



-4ttR 2 



(5.8) 



(5.9) 



which is the limiting behaviour of Eq. (5.7) (with xx' set 
equal to a value of order of the short distance cutoff, i.e. 
the lattice spacing). Note that without pinning : 

/ e i27ril0(t,O) e — i2wR<i>(Q,0) \ 

_ /gi27rR<^i,(t,0) e -i27rR^i,(0,0)\ / e -i27rR<^fl(t,0) e i27rR<^R(0,0)\ 

oc |i|- 2rf2 . 

(5.10) 

Thus pinning increases the exponent by a factor of 2. 
Perhaps surprisingly, pinning <p(0) changes correlation 
functions involving cf>. This effect ultimately arises be- 
cause 4> and 4> cannot be regarded as being independent. 
d<j)/dx is canonically conjugate to cf>. Therefore when <f> 
gets pinned we must take into account the effect on ip. 
Different results were obtained in Ref. || because this 
effect was not taken into account. In particular, a fac- 
tor of two discrepancy in the exponent occurring in the 
electron Green's function (at x — x' — 0) arose from the 
mechanism explained here. An alternative way of under- 
standing the disagreement with Ref. || can be found in 
Ref. [561 and in Ref. El. 



VI. SUMMARY 

In summary, this paper studies the effects of an impu- 
rity in a non-half-filled ID Luttinger liquid. Our numer- 
ical study for the low energy spectrum confirms previous 
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results. That is, we show that the low energy spectrum of 
a chain with an impurity flows to the spectrum of a chain 
with PBC if the interaction is attractive, and to the spec- 
trum of a chain with OBC if the interaction is repulsive. 
The behavior of a non-half-fillcd system is shown to be 
the same as that of a half-filled one and the presence of 
particle-hole symmetry is not essential for this behaviour. 
Therefore, we numerically confirm that the properties of 
a generic ID Luttinger liquid with an impurity are indeed 
those predicted by Kane and Fisher B. 
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TABLE II. Conformal field theory prediction of d j for 
PBC and OBC fixed poi n ts ( see definition in Eq. (3^5) and 
predictions in Eqs. (3.3, i-4)). 8 is the forward scattering 
phase shift. The expressions for PBC are different for a = 0.5 
and a = —0.5 . 



a 
d 
di 
d 2 
d 3 



4tv 
1 

17W 



^ 2 (1-|) 2 
t:R 2 {1 



2 
S\2 



{-\,PBC) 

irR 2 (^) 2 
7rfl 2 (l-i) 2 
1+ttR 2 (1) 2 

1+^2(1*1)2 



{±\,OBC) 

^ 2 (|) 2 
±+nR 2 (l-i) 2 

^ 2 (|) 2 
fk+^ 2 (H) 2 



TABLE III. Properties of M = model. 



a 




.5 


-.5 


L(E X - 


E ) 


W3/2 




e 2 (b = 


0) 


3/2 


3 




1) 


9/4 


G 


e 3 (b = 


o) 


5/2 


4 


e 3 (b = 


1) 


13/4 


7 



7 



TABLE IV. The scaling exponent 7 for small deviation 
from fixed points (a, fe*). 7 > and 7 < means the en- 
ergy spectrum flows away from and flows to the fixed points, 
respectively. 



(a, b-. 

1 



■) (i 1 ) 

1 - l/47ri? 2 
0.17 



(1,0) 
1 - 4vri? 2 
-0.20 



K,l) 



2(1 



0.78 



(-|,o) 

) 1 - 47ri? 2 
0.28 
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fig. 1 Qin et. al. 
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FIG. 



0.02 
1. For a = 



0.04 

0.5, b 



0.06 
1/L 

and b 



0.08 



0.1 



0.12 



1, the di defined 
in Eq.(3.5) is plotted for chain lengths L = 8, 16, . . . , 64 as 



a function of 1/L. We extrapolated the di of chain lengths 
L — 32 to 52 to L — > 00 limit by polynomials of 1/L to the 
second order. The extrapolated values for di are in agreement 
with the conformal field theory predictions given in Table. II, 
with 5 = for b = 1 and 5 = -0.29 for b = 0. 
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1.2 
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fig.2 Qin et. al. 



PBC(b=1): d 



d2 -B- 
d 3 x 



OBC (b=0): d 

d-| -iis- 
tb o- 
CI3 -t- 



o <>-o- -0 



t ! f -y 1 A A A 



0.02 0.04 

FIG. 2. For a = -0.5, 6 



0.06 
1/L 

and b 



0.08 



0.1 



0.12 



1, the di defined 
in Eq.(3.5) is plotted for chain lengths L = 8, 16, . . . , 64 as 



a function of 1/L. We extrapolated the di by data of chain 
lengths L = 32 to 52 to L — > 00 limit by polynomials of 1/L 
to the second order. The extrapolated values for di are in 
agreement with the conformal field theory predictions given 
in Table.II, with S = for b = 1 and S = -0.19 for 6 = 0. 
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fig.3 Qin et. al. 
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fig.5 Qin et. al. 



0.06 

1/L 



FIG. 3. The orthogonality exponent a(L) and ground 
state scal ing dimension do(L) for OBC chains, (see Eq.(4.7) 
and Eq. (4.8)), are extrapolated to L — > oo limit by polyno- 
mials of 1/L to the second order. We show the extrapolated 
values for a and do are equal as conformal field theory pre- 
dicts. 

fig. 4 Qin et. al. 



-0.5 



linear fitting lines: 



a=0.5, b=0.9 
a=0.5, b=0.1 
a=-0.5, b=0.9 
a=-0.5, b=0.1 
0.17-4.27/ln L 
-0.20-3.1 3/ln L 
-0.78-2.49/ln L 
0.28-3.1 5/ln L 



■**>- 



FIG. 4. The scaling exp onent 7 for excitations energy near 
fixed points, see Eq. (4.9), is given for chains with lengths 
around 50. The calculated 7 are in agreement with the theo- 
retical values listed in Table IV. 
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PBC fixed point + 
OBC fixed point □ 
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FIG. 5. The renormalization group' 'flow of the excita- 
tion energy to OBC (b — 0) fixed point for a — 0.5 is 
qualitatively shown in this figure for the range L = 20 to 
50. The line segments from left to right of the plot are for 
b = 0.05, 0.15, 0.95, respectively, and the arrow indicates 
how excitation energy A(0.5, b) flows when L increases. 



fig.6 Qin et. al. 
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FIG. 6. The renormalization group flow of the excita- 
tion energy to PBC (b = 1) fixed point for a = —0.5 is 
qualitatively shown in this figure for the range L = 20 to 
50. The line segments from left to right of the plot are for 
b = 0.05, 0.15, 0.95, respectively, and the arrow indicates 
how excitation energy A(— 0.5,6) flows when L increases. 
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fig.7 Qin et. al. 
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FIG. 7. The orthogonality exponent a{L) and ground 

state scaling dimen sion do(L) f or i mpurity chains with 

a = 0.5, (see Eqs. ( p~7j ) and Eq.(4.8)), are plotted in this 



figure. The line segments from left to right of the plot are for 
b = 0.0, 0.05, 0.15, 0.95, 1.0, respectively. The figure shows 
qualitatively how a(L) and do(L) flow to OBC fixed point as 
the chain length increases. 



fig.8 Qin et. al. 




arctanf L u ^° tan(b7l 12) ] / Jt 
FIG. 8. The orthogonality exponent a{L) and ground 

state scaling dimension do(L) for i mpu rity chains with 

a = -0.5, (see Eqs. ( ffTjj ) and Eq. (4.8)), are plotted in 

this figure. The line segments from left to right of the plot 

are for b = 0.0, 0.05, 0.15, 0.95, 1.0, respectively. The figure 

shows qualitatively how a(L) and do(L) flow to OBC fixed 

point as the chain length increases. 
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